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Abstract 

Some problems related to the structure of higher terms of the e-expansion of Feynman 
diagrams are discussed. 
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1 Introduction 



Present accuracy of experimental results |T[ demands calculations of higher-order (two-, 
three- and, sometimes, four-loop) radiative corrections. Within dimensional regulariza- 
tion [0, with the space-time dimension n = 4 — 2s, one needs to evaluate higher terms 
of the e-expansion of multiloop integrals. Dimensional regularization allows one to simulta- 
neously regulate the ultraviolet (UV) and infrared (IR) singularities. The highest power of 
1/e contributing to the UV-singular part is connected with the number of loops L. At the 
same time, IR-singularities or the reduction techniques || |[] may introduce additional pow- 
ers of 1/e. For this reason, the higher-order terms of the e-expansion of one- and two-loop 
diagrams become very important for multiloop calculations. In some cases, one can derive 
exact analytical results valid for an arbitrary space-time dimension n, usually in terms of 
various hypergeometric functions [[| |6| . However, even in these cases it may be very difficult 
to extract from these functions the analytical form of the coefficients of the e-expansion. 

In this paper, we review some recent results devoted to construction of higher-order 
coefficients of the e-expansion of (mainly) the single-scale diagrams, whose only mass scale, 
M, can be easily factorized as (M 2 )~ Le , L being the number of loops. After factorization, the 
coefficients in e for this diagram are nothing but a combination of some (in general, irrational) 
numbers. In some cases, these numbers can be numerically calculated for each order of e, 
with a very high accuracy 0. Nevertheless, only the exact analytical result would allow 
one to check important properties of quantum field theory, like gauge invariance of physical 
observables and the cancellation of the UV-poles ||. Among the physically-important cases 
belonging to this class are renormalization group (RG) calculations within MS scheme 0, 
on-shell calculations within QED and QCD JTD] and some low-energy electroweak processes, 
where all external (internal) parameters are small with respect to internal (external) ones 
TIL 



2 Algebraic classification 



Calculation of the RG functions, like /^-function or anomalous dimension 7, within MS 
scheme can be reduced to the corresponding propagator-type massless diagrams. As it was 
shown in ||, all three- loop RG functions in arbitrary renormalizable models are expressible 
only in terms of (^-functions. The four-loop results for QCD and five-loop /^-function for 
</> 4 -model confirm this theorem || . At the higher orders (six and more loops) new, non-zeta 



terms appear [[L3j. They can be expressed in terms of the multiple Euler-Zagier sums [14 



C(si,...,s fc ; 01, . . . , a k ) = IT 



Si 



ni>n 2 >...>n k >0 j=l 

where Oj = ±1 and Sj > 0. We will also use a short-hand notation 

( Sl ,...,s k = COi, • • • , s k ; 1, . . . , 1) , U atb = ((a, b; -1, -1) . 

For sums of the type (|I]), the weight j can be defined as X4L1 s i> whereas the value of k can 
be associated with the depth (see in [|TJ|, [T^]). For lower cases, these sums correspond to the 
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ordinary (^-functions. The six-loop /5-function for the scalar model contains £5,3 (this was 
recently confirmed in |17|). At the seven- loop order, a new transcendental number arises [ |i~3H , 
C3,5,3- Note that £5,3 and another constant, £7,3, appear in the calculation of anomalous di- 
mensions at 0(1/ 'N 3 ) in the large- iV limit pf. These results were confirmed and generalized 
where it was demonstrated that in odd dimension 3 
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2e the counterterms contain 
other constants U a ^. A remarkable property of all these constants is their connection with 



knots |13|, [18], |5D| . All these results confirm miraculous "link" between QFT and the knot 
theory |2(|: some Feynman diagrams can be connected with certain knots, so that the val- 
ues (Euler-Zagier sums) of Feynman diagrams are also associated with knots. A detailed 



description of the alternating (non-alternating) Euler-Zagier sums is presented in Ref. [|15 



All sums of the weight j and the depth k, satisfying condition 1 < k < j < 7, are reduced to 



12 basic elements and their products. These irreducible numbers may be taken as |T5j 



ln2, ir 2 , Cs, Cs, Cv , t/5,1 , C(5, 1,1; 1,1,-1), C(3, 3,1; 1,1,-1) , 




M - 1) ( 



7T 



(2) 



.7=4,5,6,7 



3,1 



2L14 ( ~ ) + — J2 In 4 2 + \C,2 In 2 2 . It is interesting that also arises 



In particular, U; 

in three-loop massive calculations in QED and QCD fTO , [T~2|| within the on-shell scheme 

Standard Model calculations produce diagrams with other mass distributions. Two-loop 
vacuum diagrams with equal masses 0, ^l[ yield the transcendental number CI2 ffj , where 

CL; (0) is the Clausen function ||, Chi (0) = Im Li 2Z (e ie ) and Cl 2 /+i (6) = Re Li 2m (e ie ). 
The same constant Cl 2 (j^j appears in the one-loop off-shell three-point diagram with mass- 
less internal lines, in the symmetric case when all external momenta squared are equal [^3| . 
A generalization to the L-loop ladder case yields C1 2 l (|J p4 . 



At the three- loop level, however, new constants appear |fL2j| . For their classification, 
Broadhurst [PJ has introduced the "sixth root of unity" basis connected with 
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3=1 n i 



(3) 



where A = e 17r ' 3 and pj G {0, 1, 2, 3, 4, 5}. For pj e {0, 3} it coincides with the Euler-Zagier 
sums ([!]). For the lowest weights j, the bases consist of the following elements: 



1 

2 



7T, In 2, In 3. 

7T 2 , In 2 2, In 2 3, Cl 2 (f) , Li 2 (±) , In 2 In 3, rrln2, rrln3. 
3 : 7T Cl 2 (f ) , 7T 2 In 2, tt 2 In 3, Cl 2 (f ) In 2, Cl 2 (f ) In 3, rr 3 , Li 2 (|) In 2, Li 2 (j) In 3, In 3 2, 
In 2 2 In 3, In 2 In 2 3, In 3 3, vrLi 2 (±) , vrln 2 2, vrln 2 3, vr In 2 In 3, ( 3 , Li 3 (^) , Li 3 (f) . 

Unfortunately, a large number of elements (more than 4000) makes it difficult to define the 
complete basis of "sixth root of unity" at the weight 4. So far, only the cases with the depth 
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k < 2 have been examined []16 |. One of the remarkable results of Ref. Jl6[ is that all finite 
parts of three-loop vacuum integrals without sub divergences, with an arbitrary distribution 
of massive and massless lines, can be expressed in terms of four weight-4 constants: (4, 



U 3; x, and V 3t i 
2 



Z^p>fc>0 „3k 



cos 



Cl 2 (f 
Note that 

internal lines are massless, whereas all external momenta squared are off shell and equal. 



CI2 (f ) 



(fvrA;) , which is an essentially new constant, 
also appears in the two-loop non-planar three-point diagram |25j , when 



In Refs. [pq, |27| , the central binomial sums were considered, S'(a) = J2^= 



(n!) 
(2n)] 



In 



particular, in [26] one-fold integral representation was established. It was shown that these 
sums are connected with the multi-dimensional polylogarithm E3] of "sixth root of unity" , 



Liai,...,a fc (A) 
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a k 1 
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and Euler-Zagier sums. 



3 Geometrical approach 

To predict types of functions (and the values of their arguments) which may appear in higher 
orders of the e-expansion, a geometrical approach |29| happens to be very useful. Using this 
approach, the results for all terms of the e-expansion have been obtained for the one-loop two- 
point function with arbitrary masses 0, . Moreover, all terms have been also obtained for 
the ^-expansion of one-loop three-point integrals with massless internal lines and arbitrary 
(off-shell) external momenta and two- loop vacuum diagrams with arbitrary masses [Rj], RJJ, 
which are related to each other, due to the magic connection J33|. All these results have 



been represented in terms of the log-sine integrals (see in |22j and below), whose angular 
arguments have a rather transparent geometrical interpretation (angles of certain triangles). 
For instance, for the one-loop two-point function with masses mi and 777.2 the relevant angles 

are (see in |2J, |3TJ) r 01 = arccos ( 2mi ]/^ 2 ) and r 02 = arccos { 2 m 2 VW )' where k is the 



external momentum. For single-scale diagram with mi = 7712 = 777 and k 2 = m 2 , each of 
these angles is equal to |. In more complicated cases, like, e.g., the three-point function 
with general values of the momenta and masses, an arbitrary term of the ^-expansion can 
be represented in terms of one-fold angular integrals whose parameters can be related to the 
angles accociated with a three-dimensional simplex pi RQ]. 



In Ref. fl34|, the on-shell values of two- loop massive propagator-type integrals have been 
studied, and it was observed that the finite (as e — > 0) parts of all such integrals without 
sub divergences can be expressed in terms of three weight-3 constants, two for the real part, 
(3 and 7rCl2 (f), and one for the imaginary part, tt( 2 - 

It is natural to assume that the elements occurring in the coefficients of the e-expansion 
are connected with some properties of diagrams, like the structure of the cut(s). In this case, 
the numbers of elements for higher-order weights may be essentially reduced. Basing on this 
conjecture, in Ref. an ansatz was elaborated for constructing the "irrationalities" occur- 
ring in the ^-expansion of single-scale diagrams involving cut(s) with two massive particles. 
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This construction is closely related to the geometrically-inspired all-order e-expansion of the 
one-loop propagator-type diagrams 13 H, which was also used to fix the normalization 
factor The procedure of constructing the ansatz is as follows: for each given weight j 
the set {bj} of the basic transcendental numbers contains (i) all products of the lower- weight 
elements k = 1, 2, ■ ■ • , j — 1 and (ii) a set of new (non-factorizable) elements {bj}, 

which are associated with the quantities arising in the real and imaginary parts of the poly- 
logarithms Lij (V 6 ^ and Lij ^1 — e l6 *^, with 9 = | or 9 = They can be expressed in terms 
of the Clausen function Clj (9), log-sine integrals Lsj (9) and generalized log-sine integrals 
Lsf } (9), defined as (see in O) 



Lsf ) (9) 



ln J 



-fc-i 



2sin§ 



L Sj (9) = Lsf (9) . 



(5) 



Note that Ls2 {9) = CI2 (9). Therefore, in our case the non-factorizable part of the basis can 
be expressed in terms of the generalized log-sine integrals of two possible angles, 9 = | and 
9 = To establish the basis of the weight 4 within such an ansatz, one needs to analyze 
no more than 100 elements. It should be noted that this basis is not uniquely defined, since 
there are several relations between polylogarithmic functions Clj (9), Lsj (9) and Ls^ (9) 
of these arguments. After excluding all linearly-dependent terms (see Appendix A of fl36f), 
the basis contains the following non-factorizable constants: Lsj li=3 4 5' -^ s i (f)lj=2 4 5' 

(t) 1.7=4,5 an d (x) - This se ^ °^ e l emen t s is called the odd basis [05 |. The numerical 
values of these constants are given in Appendix A of [35|j. In particular, it was found 
that the constant V3 1 can be expressed in terms of the weight-4 elements of the odd basis, 



V, 



3,1 



Cl 2 (f)] -ivrLsa 



2tt\ , 13 r ln o 259/- , 3t „(1) ( 2tt\ 

By analogy with the odd basis introduced in ||35|| , it is possible to consider the even basis 
connected with the angles | and it. Apart from the well-known elements ir, ln2, Q and 

the Catalan's constant G = CI2 (|) = LS2 (f), this basis also contains (up to the weight 5) 

Lij (~) \- =4 5 (see also in [0), Lsj (-} | J=3 4 5 , CI4 f |) and Ls^ (j^J . An example of a physical 

calculation where the constant LS3 (|) arises is given in []37|]. Instead of Li 4 (J^J and Li 5 Qj, 

one could take, e.g., Ls 4 ^ f|) and LS5 f|) (or Ls 4 ^ (7r) and Ls^ (ir)) (see Appendix A of 
|36|| ). The constructed odd and even bases have an interesting property: the number Nj of 
the basic irrational constants of a weight j satisfies a simple "empirical" relation Nj = 2 3 ', 
which has been checked up to weight 4. The situation with weight-5 bases is discussed 
below. 



4 Searching for new elements 



A natural generalization of the ansatz proposed in []35] could be the inclusion of the gener- 
alized (Nielsen) polylogarithms f38| , S aj b(z), where again we consider the cases z = e ld and 
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z = 1 — e l9 , with 9 = | or It is easy to show that 



i a (-l^ 



d0 (6 - 4>) 



a-1 



hi 



2sin§ 



+ S a , b (1) 



(a - 1)!6! Jo 

reduces to a combination of the generalized log-sine integrals (|5|). The function S a> b (l — e 10 ) 
can be reduced to S a ^ (e l6 ) plus products of lower-order S a ^ (e l9 ) and powers of In ^2 sin 



Therefore, such Nielsen polylogarithms do not generate new elements. Note that S a ,b (1) is 
a combination of (^-functions |38| . 

Another possibility is the multi-dimensional polylogarithm (0). Performing PSLQ anal- 
ysis with an accuracy of 100 decimals we established that, up to weight 5, all constants 
([§) are expressible in terms of the odd basis. 

Let us consider the following hypergeometric function: 



p+i-fp 



1 + a x e, . . . , 1 + a P+1 e 
+ be, 2 + cie, . . . , 2 + c R e, 1 + die, . . . , 1 + d P ^ R _iS 



(6) 



where < z < 4. The one-, two- and three- loop Feynman diagrams which are expressible 



via hypergeometric function of this type are given in 



The e-expansion of the function 



can be written as (for details, see Appendix B of Ref. Pq|) 



2(i + 2te)r* 

11(1 + Cie) 



.i=i 



E 



(j!) 2 J 



R+l 



1-6 



5iTi + — + 2bS 1 
J 



T2 O + 



£>1 



+2bS 1 S 1 T 1 + — S-jTi + 26S' 1 + 2b l [S 2 + Sf) + -S 2 T 2 + -SfT{ 



1 



2j 
+ 0(e 3 



(7) 



where A k = Ef^ 1 a*, C k = £f =1 £ fe = E^" 1 and T fe = C fc + D fe - A k - 6 fc . It is 
implied that C k = and Ilili(l+ C i£) — 1 when R = 0, and .D^ = when P — R—l = 0. Here 
and below, we also use the short-hand notations S a = S^n — 1) and S a = S a {2n — 1), where 
S a (n) = Ej=i J _a is the harmonic sum. In this way, the e-expansion of the hypergeometric 
function @ is reduced to series of the following type: 

w = E SS^)" • • • ( s J ip ■ ■ ■ ( 8 ) 



J ai ,...,a p ; bi,...,b q ;k\ 



£i(27i)!n* 



For p = g = 0, the analytical result is available [[10 

fc— 2 



00 (n!) 2 ^ n 



E 



E 



-2V 



-(lnz) fc - 2 ^Lsg 2 



9 Z = 2 arcsin 



^(2n)!n* (k - 2 - " 2 

For the p + g = 1 terms, the following one-fold integral representation can be constructed: 

2 



(9) 



oo / |\2 n 

E|^!^ [ ^ l(n - 1) + 5l(2r2 - 1)] 

x {±Ls 2 (9) + ±01n (2 sin |) - (1 + a) 



(k-2) 
Ls 2 Ok + < 



jf " d9 In z - 2 In (2 sin |) 
) + #ln(2cos§)]} , 



fc-2 



(10) 
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where k>2,0<z<4 and a is an arbitrary constant. For all other types of the sums 



the results are available only for particular values z = 1,2,3 H36| , (40f . It should be noted 
that not all sums @ can be expressed separately in terms of the even and odd bases, and 
vice versa, not all basis elements are expressible in terms of the binomial sums. 

Testing the e-expansion of a one-loop three-point diagram contributing to the H — > 77 



decay (see details in [|3~6jl ) we revealed that new elements X5 an d Xs appear at the weight-5 
level in the odd and even bases, respectively. They can be presented as special types of the 
inverse binomial sum (|8T), 



00 (n\) 2 1 



X5 



00 (n\) 2 2 



X5 



'ID 



The r.h.s. of Eq. ([TI]) gives us an idea about a possible integral representation for new 
functions arising in higher orders of the e-expansion. Apart from Ls^ (#), we meet 



Lscjj (6) 



ln J 



2sin§ 



ln J " 



2cosf 



(12) 



Some properties of these functions were studied in Appendix A of |BB[ . In particular, a careful 
PSLQ analysis showed that a new element for the even basis should be added, Lsc^^ ). In 
this way we restore, accidentally, the "empirical" relation Nj = 2 J for the j = 5 level of the 
odd and even bases. However, the class of possible new functions is not restricted only by 
Lscjj (8) (|T2|). We have not tested the harmonic polylogarithms |^1[ and their multiple- value 
generalization |42| , which may be related to the e-expansion of certain Appell functions (see, 



e.g., in [@). Another possibility is to test all types of the integrals occurring in fllPl): 



d0 fc In 1 



2 sin ■ 



2 cos ■ 



dcj) Ls 2 ((p) W 



2 sin 



2 cos • 



(13) 



All these issues require a more careful analysis. We present here one-fold integrals related to 
our elements X5 an d X5 obtained by using a PSLQ analysis with an accuracy of 200 decimals: 



2vr/3 



Ls 2 (0) In 2 (2 sin f ) = -§-^2^2 (f ) + HC2C3 + HCs + H 7 ^ (f 
-^-Ls 4 (f ) - !Ls 2 (f) Ls 3 (f ) + |Ls« (f ) - ± X5 



0.744148409838194515377332007. . . 



(14) 



fj 2 d<P Ls 2 (0) In 2 (2 sin f ) = -i In 5 2 + ^(2 In 3 2 - ^Cs In 2 2 - §Li 4 (|) In 2 - §Li 5 (| 



475 



48 1 48 sz 128 Sd 8 4 V27 8 °V2 

i024-C 2 C3 + HC5- ^vrC 2 Ls 2 (f ) - ±nLs, (f ) + ^vrCl 4 (f ) - |Ls 2 (f ) Ls 3 (f 
: 0.656546903118158978863954386 ... (15) 



48 X5 



5 Concluding remarks 

Let us return to the idea about the connection between the structures occurring in the 
£-expansion and the structure of the cut(s). Careful analysis of existing results [13] for 
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the single-mass-scale diagrams allows us to formulate the following conjecture jMj: the e- 
expansion of diagrams with cuts involving 0, 1 or 3 massive lines (or their combinations) are 
expressible in terms of the even basis, whereas the diagrams with a cut involving 2 massive 
lines (or its mixing with a or 1 massive cut) produce elements of the odd basis. For all other 
cases the problem is open. Recent results [fE| confirm this conjecture. We would expect that 
the four-loop numerical results in QED j46| can be represented in terms of the even basis. 

Another interesting example is the one-loop triangle diagram J3 with three massive lines 
with rrii = m, when all scales (masses and momenta) are equal to each other. The geometrical 
approach [TM BO] yields the following representation: 



1,1,,. 



1: m) 



17T 



2-e 



m 



2+2e 



r(i+ £ )2 1 - £ 3 1/2H 



3=0 



tt/3 



d01n i+1 l 14 



■ cos^ 



(16) 



where (see also in Refs. pTL H3, 



tt/3 



In 1 



8 cos 2 



7T . 



f 1 + ^;) 

Jo \ 8 cos <p J 



In2-Cl 2 (f) 
|Cl 2 (f)ln2 + |vrln 2 2 + |7rLi 2 (i 



(17) 



|Ls 3 (jf ) 



+2Im 



Li 3 {\t 



nr/3 



54^2 • 



One can see that the e-term is expressed in terms of "sixth root of unity" basis |16| . However, 
it is not sufficient to have only the elements from the odd/ even bases. To explain this 
situation, we should recall that this vertex diagram also possesses an anomalous threshold. 
If we put all p 2 = p 2 then this threshold would correspond to p 2 /m 2 = 3. In other words, 
this triangle diagram has a mixture of an anomalous cut corresponding to p 2 /m 2 = 3 and 
two-particle cuts (with respect to each leg) at p 2 /m 2 = 4. A similar situation is also observed 
for the case when p 2 /m 2 = 4 and p 2 /m 2 = 9 cuts are present |49 . 



Finally, let us summarize the situation with the even and odd bases. Although our con- 
struction is incomplete, it has allowed us to obtain new results for several one-, two- and 
three-loop master integrals implemented in different packages ||50|| . In this way, we found 



several interesting relations between generalized log-sine integrals and multiple zeta func- 
tions. Moreover, a new relation between 3F2 and 2-F1 hypergeometric functions of argument 



was established 136] . 
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